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Abstract 

The shock jump conditions for the Euler equations in their primitive form are 
derived by using generalized functions. The shock profiles for specific volume, speed, 
and pressure are shown to be the same, however density has a different shock profile. 
Careful study of the equations that govern the entropy shows that the inviscid entropy 
profile has a local maximum within the shock layer. We demonstrate that because of 
this phenomenon, the entropy propagation equation cannot be used as a conservation 
law. 
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1 Introduction 


A consequence of the nonlinearity of the equations of motion is the steeping of compression 
waves into a shock wave. Within the shock layer, the gradients of velocity and tempera- 
ture become large, and irreversible thermodynamic processes caused by friction and heat 
conduction become dominant. At high Reynolds numbers, the shock- layer thickness is of 
the order of several mean free paths; for all practical purposes, the shock layer can be 
represented as a mathematical abstraction that corresponds to a surface across which the 
flow variables experience a sudden jump. Away from this discontinuity surface, viscous and 
heat conduction effects are usually negligible and the inviscid equations of motion model 
the flow well. Remarkably the information needed to account for the final outcome of the 
irreversible processes that take place within the shock layer is contained in the inviscid 
equations. 

The study of shock waves is 150 years old. The jump conditions satisfied by the conserva- 
tion of mass and momentum were discovered by Stokes [5] in the middle of the 19th century. 
Stokes’ excitement at making this discovery is evident in his paper: “These conclusions cer- 
tainly seem sufficiently startling; yet a still more extraordinary result. ..the result, however, 
is so strange...” The shock jump condition associated with the conservation of energy was 
implicit in an investigation conducted by Rankine [7] in 1870; however a precise exposition 
was not made until the work of Hugoniot [5] in 1889. The increase in entropy across a shock 
was a more difficult concept to grasp. The leading fluid dynamicists in England (Stokes, 
Kelvin, and Rayleigh) questioned the validity of the shock discontinuity because it violated 
the conservation of entropy. The correct principles were not well understood, and did not 
appear in their present form until around 1915. 

Detailed studies of the viscous shock layer emerged several years later with the work of 
Becker [1], who solved exactly the one-dimensional equations of a real fluid. In a related 
study, Morduchow and Libby [6] found the exact entropy distribution across the shock 
layer of a viscous heat-conducting gas. Morduchow and Libby observed, see Figure 1, that 
the entropy, unlike the other flow variables that behave monotonically, increases through 
the shock layer until it reaches a maximum at the center of the layer and then decreases 
to its expected value on the other side of the shock. Morduchow and Libby explained 
this phenomenon as follows: “It may at first sight appear that the recovery of mechanical 
energy on the downstream side of the center of the wave thus indicated by this solution 
would violate the second law of thermodynamics... However, the second law applies to an 
entire system-that is, to the end points-and permits energy recovery in separate sections 
thereof. The negative entropy here might also be interpreted as physical effects that are 
not taken into account by the governing equations...” 

Today, the shock jump conditions are obtained for the inviscid equations by casting them 
in their integral conservation form. A brief derivation, based on this standard procedure, is 
given section 2. However, the purpose of this work is to show that the shock jump conditions 
can be derived from the primitive differential form of the equations. The genesis of the 
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of the analysis presented here was contained in an unpublished work of Gino Moretti written 
in the early 1970’s. The significance of this work is primarily that it demonstrates how 
to obtain the shock jump conditions for equations that cannot be cast in a conservation 
integral. Similar work has been presented by Colombeau in reference 2. An interesting 
consequence of this exposition is a better understanding of why the entropy equation does 
not yield the proper jump. 


2 Standard Shock Jump Analysis 


The derivation of the jump conditions across a shock associated with the Euler equations 
is well known. See, for example reference 9. The derivation is included here so that it can 
be contrasted with the “nonstandard” approach introduced in section 3. 

The one-dimensional conservation laws for the inviscid flow of a perfect gas are 


J [pt + (pu) x 


dx — 0 


pu)t + (pu 2 + p) x 


dx = 0 


J \{pE) t + (pu H) x ) dx = 0 


( 1 ) 


where p is the density, p is the pressure, E is the specific total energy, H is the specific 
total enthalpy such that 

H = E + - 
P 

and u is the velocity of the gas. In the standard analysis for the shock jump conditions, we 
weaken the usual smoothness requirements associated with the classical notion of a function 
by introducing the concept of a weak or generalized solution. Basically, the integrand of 
equation 1 is multiplied by a test function that is at least C 1 smooth and has compact 
support. Then, w r e integrate over space and time in the neighborhood of the shock and use 
integration by parts to move differentiation from the discontinuous fluid variables onto the 
smooth test function. Thus let equation 1 be symbolically represented by 


/ ( u t + F x ) 


dx — 0 


( 2 ) 


where 

and 


U = (p, pu, p E) 1 


F = (pu, pu 2 -f p, pu H) t 
and let the initial conditions be given by 


U(x,0) = U o (x) 
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Let <p be a test function that is continuously differentiable and has compact support. Con- 
sider the domain D around the shock E defined in the rectangle 0 < t < t\ and a < x < b, 
see figure 2. Let <f) be zero outside of D and on its boundary. Multiply the integrand of 
equation 2 by <f>, integrate over x and t, and use integration by parts to obtain 

/ f {U 6 t + F d> x ) dx dt — 0 (3) 

Jd Jt > o 

We say that U is a weak solution of the initial value problem 

Vt + F x = 0 

with initial data Uq if equation 3 holds for all differentiable test functions <p with compact 
support. 


Let Di be the subset of D on the left of E and D r be the subset on the right of E as in 
figure 2. Assume that U is differentiable everywhere except across E; hence, on Di with 
the divergence theorem, we find 

[ [ (U <k + F <f> x ) dx dt = / [ [(U<f>)t + (F(j)) x } dxdt= [ d> i-Udx + Fdt) (4) 

JDi Jt> 0 Jd { Jt > 0 JdD 

and similarly for D r . Because <p is zero on the boundary of Z), the line integrals are only 
nonzero along the shock E. Let the shock be defined by x(t) and Ui be the value of U on 
the left of the shock; similarly, let U r be the value of U on the right side of the shock. Then 
by using equation 4 and the equivalent expression on the right of E, we obtain 

J 6{-{U\dx + [F]dt) =0 

where (//] = U r — Ut and [F] = F{U r ) — F{Ui). Because o is arbitrary, 

c| U] = [F] (5) 

along S. where c = dx/dt is the speed of the shock. Equation (5) results in the following 
Rankine-Hugoniot (R-H) jump conditions: 

\p(c-u)} = 0 

[pu{c- u)] + [p] = 0 (6) 

[pE(c - u)] + [pu] = 0 

One solution to equation 6 corresponds to no mass flow across the discontinuity and leads 
to the conditions across a slip line. The other solution results in jumps in pressure, density, 
and velocity. After some manipulation, the R-H jumps can be expressed as 


pi Ui — Pr fir 
Pi + {pu 2 )i = 


Pr + [pU )r 


1 2 , l-2_ 1 

+ 2 Ul ~ 


«r + 


(7) 
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where u = u — c ? a is the speed of sound, and 7 is the ratio of specific heats. The above 
relations indicate that 

Pr ~ Pi 


n r ui — 

Pr - Pi 

which is known as Prandtl’s relation, and 

Ul_ _ (7 + l)Pr + (7 - 1 )pi 
U r (7 - l)p T + (7 + l)pi 

These results imply that the entropy jumps across a shock. Its jump is given by 


[S\ = In v ~ 

Pr 


In ^ 

Pi 


( 8 ) 


(9) 


( 10 ) 


Although the entropy propagation equation can be expressed in the form of a conservation 
law as 

' (pS), + (puS%]dx = 0 (11) 


it cannot be used to obtain the correct entropy jump across a shock wave. 


3 Nonstandard Shock Jump Analysis 

The shock jump conditions can be derived without relying on the integral conservation 
laws. The importance of this method is threefold. First, this method provides a means for 
determining the jump conditions for physical laws that cannot be expressed in conservation 
form. Second, it will ultimately lead to an understanding of the nature of the entropy 
structure across an inviscid shock. Thirdly, it may suggest how to derive shock-capturing 
algorithms with proper jumps from the nonconservation form of the Euler equations. 
Consider the following system of equations: 

V, — vu x + uv x = 0 

u t + uu x + vp x = 0 
Pt + up x + 7 pu x = 0 

where v = 1 j p. The reason for using v instead of p will become clear in section 4. 


We look for solutions to v, u. and p of the form 


V = if, + [v] H(Z) 

(15) 

U - u; + [u]/t(f) 

(16) 

P = Pi + [p}L{0 

(17) 

here for now we only require that 


TI T , T ( 0 for x — > -oo 

H , A , L = < , 

( 1 lor x — > oc 

(18) 


( 12 ) 

(13) 

(14) 
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£ = x — ct, and [ w ] = w r — w t . The functions H , K, and L provide a description of the 
shock profile or structure with end conditions wi at x = — oo and uy at x = oo, where w 
stands for e, u , and p. 

Consider equation 12 and introduce equations 15 and 16: 


— c[v]H' — (vi + [■ v]H)[u]K ' + ( ui + [u]/\)[u]/T = 0 
We can then rewrite equation 19 as 


dH H 


Vl 


dl\ a + I\ [e](a + K) 


where 


By integrating equation 20 we obtain 


a = 


Vl 


Ul — c 

"ET 


H — t— r + b( a + A ) 

[v] 

where b is a constant of integration. Now, for x — ► — oo, both H and K 
the constant of integration is 


b = 


t>]a 


and 


» = A- 

[u] a 


In the same way as for x — > oo, both H and K — » 1,; therefore, 

vi 


[u]a 


= 1 


This last, relation, together with equation 21, gives (with u = u — c ) 


Ur _ pi_ 
Ul Pr 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 

0,; therefore, 

(23) 

(24) 

(25) 


(26) 


which is the R-H jump for the conservation of mass equation found in section 2 (equation 
6). By using equations 24 and 25, we find that 


H = K 

Now consider equation 13; if we introduce equations 15 - 17, then 

- c[u]I\' + ( ui + [u]A')[u]/C ' + (vi + [v]H)[p]L' = 0 


(27' 


(28) 
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With equations 27 and 25, we can rewrite equation 28 as 


We can integrate 


dL [u ] 2 

dK + [v]\p] 


L + ±±-K = d 

Mb] 


where d is a constant of integration. As K — ► 0, L 
Hence, for A' — » 1, L — ► 1, 

bl 2 


Mb] 


i 


This reduces to Prandtl’s relation: 


U r U\ 


Pr~ Pi 


Pt Pi 

If we use equation 31 in equation 30, we get 

K = L 


( 29 ) 


(30) 


0, we can conclude that d — 0. 


(31) 


(32) 


(33) 


From the first two of equations 12 - 14, we find that the functions H , K, and L must be 
the same to obtain solutions as in equations 15 - 17. 

Now consider equation 14 and introduce equations 16 and 17: 

\p\( ui - c)L ' + [u]K f + [u][p](KL f + 1 LK f ) = 0 (34) 

However, we see from equation 33 that K = L such that 

b](«/ - C )K' + 'fpi[u]I< > + [u]b] (7 + V)KK' = 0 (35) 


Let us integrate equation 35 from x = — oc to r = 00 as 

{bl(«/ - c ) + 7P/M) / dK + b]b] (7 + 1) / KdK = 0 (36) 

Jo Jo 

which can be reduced to the final R-H jump condition: 

Ul_ _ (7 + j >r + (7 - ])pi , 37 x 

Ur (7 - l)Pr + (7 + l)pi 

This analysis has shown that if H = K = L (i.e., if the shock profiles for v,u, and p 
are identical), then the full set of correct R-H jump conditions can be recovered from the 
equations in primitive form. 
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4 Multiplication of Discontinuous Solutions 


Consider equation 35 rewritten as 

{[p](u; - c) + 7P/[tt] + [u][p](7 + 1)A'} K' = 0 


( 38 ) 


If I\ belongs to the set of C 1 functions, then the equation above only admits K = constant 
as a solution. If K is allowed to be a discontinuous solution, then K' cannot be factorized 
from equation 38. Here, we follow the mathematical construction of generalized functions 
proposed by Colombeau in reference 2 and Colombeau and Le Roux in reference 3 . The 
main advantage in using this construction is that most of the operations admissible with 
smooth functions can be defined for discontinuous functions, including differentiation. (See 
appendix A.) We restrict our attention to the Heaviside function and its derivative, the 
Dirac delta function. 

The Heaviside function is such that. 


H(x) 


0 for x < 0 

1 for x > 0 


(39) 


The Dirac delta function <5(.r) is 0 in [— oo, 0[ U ]0, oo] and is such that 

/ 8(x)dx = 1 

J -OC 

Let CqGO) denote the set of all C 00 functions on D with compact support. Given G'i(z), 
G 2 (x), and the test function ty(x) € C^°(D), if 

/ [Gj(z) - G 2 (x)} ty(x)dx = 0 

Jn 

for all 'P(z), then we say that G'i(z) and G 2 {x) are associated and write Gi(z) ~ G 2 {x). 
According to the definition above, G\ ~ G 2 does not imply that G 2 G\ ~ G3 G 2 , where G3 
is some other function. Consider, for example, H n (the nth power of the Heaviside function 
H). We can show that H n ~ H. but H n H' is not associated with H H' . In fact, we have 

H n ~ l H' ~ -H' (40) 

n 

Note that if we replace the associative symbol ~ in equation 40 with the equal sign, then 
we obtain, for example, H H f = |i/' and multiplication by H yields H 2 H = | H H f . Now 
by substituting an equal sign into the equation we get | H r = \\H\ which is absurd. 


In conclusion, if we replace the equal sign in equation 35 with the symbol ~ then the 
subsequent integration is fully justified. 


In another example, we present a case in which the Heaviside functions describing the shock 
are not equal because of their behavior at 0. Consider the mass conservation equation 

pt + pu x + upx — 0 (41) 
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If we seek a solution for p of the form 


P = Pi + [/>]/(£) 


then by using the same procedure above one obtains 


/ = 


a + 1 

77 A 

CL I\ 


(42) 


(43) 


where a is given by equation 21. Equation 43 is significant because it shows that the 
microscopic behavior of the Heaviside function that describes the p jump is not the same 
as that for p or u. 


5 Entropy Structure in a Shock Wave 

As pointed out in section 2, although a “conservation law” can be written for entropy, this 
law does not lead to the correct jump. Although this fact is well known, the reason why 
is not well understood. In this section, we show that the shock profile that corresponds to 
the entropy cannot be represented by a Heaviside function; hence, the entropy propagation 
equation does not yield the correct jumps. 


Consider equations 28 and 34. Multiply equation 28 by 7 and then divide by v ; divide 
equation 34 bv p and add the two equations. We obtain the following by simplifying: 


7 (Pi + \p]L)[v]H' + {vi + [v\H)\p\L' = 0 
Because H = L we have H'L ~ H L' ~ \H' , and we get the R-H jump 

Vr fail ~ l)+pr(l + 1 )] = [vi(p T ( 7 - 1 ) +Pl(l + 1 )] 
This shows that the equation 


Pt v t / p x Tr 

— + 7- + u — + 7— 

V P v 


= 0 


p v 

has a valid jump. However equation 46 may be rewritten as 

St T u S x = 0 

because dS — dp/p -f 7 dv/v. Now we look for a solution for S of the form 

S = S, + [S]T(0 


(44) 

(45) 

(46) 

(47) 

(48) 


where T is a Heaviside function. If we substitute equations 48 and 16 into equation 47, 
then we get 

{-cT + (u, + M K ) T'} [5] = 0 (49) 
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either the expression within the braces must be 0 or the jump [5] must be 0. In general, 
because no relation exists between T and A', the expression within the braces is not 0, 
hence, -we conclude that equation 49 gives the wrong jump, namely [5] = 0. In actuality, 
the problem lies elsewhere. In going from equation 46 to 47, we have gone from an equation 
with two jumps [p] and [u] to an equation with a single jump [5]. By combining the two 
equations, we have lost some information; furthermore the assumption that the solution 
can be expressed as in equation 48 is incorrect, which will be shown below. 


Consider the following. Without loss of generality, let vi = I and pi = 1 and take Si as the 
reference state for entropy. Because S by definition is 


S = lnp + 7 In v 


( 50 ) 


we obtain, given equations 15 and 17 and the fact that L = //, 


dS_ _ dS_ H , _ (1 + [v\HT 1 {(7 + l)[p]M# + [p] + 

d( dH pi’ 7 


( 51 ) 


Because the jump [p] is 


[p] = 


27 [i 


(7 + l)[v] + 2 

we can obtain the following by substituting and simplifying: 


dS_ _ {1 + 1 n orn 7 ( 7 + tji 

d£ pi’ 7 (7 + l)[i’] + 2 


( 52 ) 


( 53 ) 


This equation shows that dS/d£ ~ 0 because (1 - 2 H)H' ~ 0 see equation 40; hence, 5 
has a maximum at the origin, which is confirmed by a study of the second derivative. As a 
result, we conclude that 5 cannot be described by a simple Heaviside function but rather 
by the sum of two Heaviside functions 


S = Si + (5* - S t )T(0 + (S T ~ S*)N(0 


( 54 ) 


where S‘ is the value of 5 at the maximum. (See figure 3.) The figure shows clear that by 
using a. single equation such as equation 47 for the entropy, we cannot determine the two 
jumps that actually represent the structure of the entropy at the shock. A more technical 
proof is given in appendix B. 


6 Conclusions 

The shock jump conditions for the Euler equations in thier primitive form were derived 
using generalized functions. It was shown that the structure of the shock profile is not 
the same for all variables. A study of the entropy propagation equation showed that if 
the shock structure for the entropy is represesnted by a single Heaviside function, then the 
wrong entropy jump is obtained. It was then shown that the proper representation of the 
entropy profile requires two Heaviside functions, but not all the information required to 
specify this profile can be obtained from the entropy propagation equation. 
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Appendix 

A Some Definitions About Generalized Functions 

Here we briefly discuss generalized functions. For more details on this subject, see references 
8 and 4. 

A test function f(x) exists such that 

1 . <f{x) is C°°. 

2 . <b{x) has a compact support (i.e. , <f{x) vanishes outside of some compact interval 
[a, 6 ]). 

Furthermore, a sequence <j> n (x) of test functions converges to 0 if 

1 . for each k. the sequence of the kth derivative <j>[ k \x), 4>[ k) (x),... converges uniformly 
to 0 . 

2 . every </>„(.t) vanishes outside a given interval [a, 6 ]. 

A generalized function T is a mapping from the set V of all test functions into the real or 
complex numbers such that, if < •, • > is the internal product operator, then we have 

1 . < T,af(x) -f btl'(x) >= a < T,<j)(x) > +b < T,f.’(x) > 

2. if d> n (x) converges to 0 in the manner defined above, then < T,(f> n (x) > converges to 
0 . 

Generalized functions are useful because their derivatives are always well defined. In fact, 
we have the following definition 

< T\6 >= - < T,4> > 


B Entropy Structure Proof 

Consider equation 50 substitute equations 12 and 14. Because H = L, we have 

S = ln(l + [p}H) + 7 In (1 + [v]H) 

The equation above can be written as 

S = In (p r ) F + 7 In (v r ) G 
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where F and G are two Heaviside functions defined as 


In (1 4- \p]H) 
In (p r ) 


ln(l + [u]f/) 
In (v T ) 


We show now that F ~ G r , although they are not identical in the sense that F H' is not 
associated with GH\ We take a test function <f>(x) (defined in appendix A) and compute 
the following integrals: 

r F<f>{x)dx= r ]n ( i+ \pw < t > ( x )dx= r<f>(x)dx 

7-oc 7-oo m(p r ) 7o 

r G<t>{x)dx= r ^( x ) dx = r<t>(x)dx 

7—co 7—oo in (u r ) Jo 

from which we conclude that F ~ G in accordance with the definition of association. To 
verify whether F H 1 ~ GH\ we have 



FH’dx 


GH’dx 


[ l In (1 + \p]H) = 1 _ 1 

A) In (p r ) [p] In (p r ) 

/■ + = 1 + T -j-fi-r 

Jo In (tv) [uj in (tv) 


Thus, we conclude that F H‘ and G H 1 are not associated. Therefore, F and G must be 
considered as two locally different Heaviside functions, although they are associated. 
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3: Entropy structure across inviscid shock 


wave. 
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